Anomalous dynamics in two— and three- dimensional Heisenberg-Mattis spin glasses 
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We investigate the spectral and localization properties of unmagnetized Heisenberg-Mattis spin 
glasses, in space dimensionalities d — 2 and 3, at T = 0. We use numerical transfer-matrix methods 
combined with finite-size scaling to calculate Lyapunov exponents, and eigenvalue-counting theo- 
rems, coupled with Gaussian elimination algorithms, to evaluate densities of states. In d = 2 we find 
that all states are localized, with the localization length diverging as as energy lo ^ 0. Loga- 
rithmic corrections to density of states behave in accordance with theoretical predictions. In d = 3 
the density-of-states dependence on energy is the same as for spin waves in pure antiferromagnets, 
again in agreement with theoretical predictions, though the corresponding amplitudes differ. 
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I. INTRODUCTION 

The study of low-lying magnetic excitations in 
quenched disordered systems presents a number of chal- 
lenges. While the absence of translational invariance is 
a complicator arising in all aspects both of static and 
dynamic behavior of inhomogeneous magnets, investiga- 
tion of spin waves is made even harder because, in many 
cases of interest, the exact ground state configuration is 
not known. 

One way around the latter obstacle has been to resort 
to simplified model systems for which the exact ground 
state is known, but which nevertheless still display non- 
trivial dynamical features. Such features, it is expected, 
may shed light on the behavior of their experimentally- 
realized, rather more complex, counterparts. 

Here we deal with vector spin glasses, i. e., Heisen- 
berg spins with competing ferro- and antiferromagnetic 
interactions. It is known that the simplest realization 
of the Edwards- Anderson picture, where one has equal 
concentrations of positive and negative nearest-neighbor 
bonds of equal strength, leads (in lattices of space di- 
mensionality c? > 1) to frustration and, consequently, to 
a macroscopically degenerate (classical) ground state. 

The drawback just described does not arise in Mattis 
spin-glasses, where the Mattis transformation^ "gauges 
away" disorder effects, as far as most static aspects are 
concerned. It is known that the Mattis transformation 
does not remove the disorder effects in the dynamics of 
these so-called Heisenberg-Mattis spin glasses, which is 
non-trivial. Indeed, investigations of spin-wave propa- 
gation in such system a^ii^i^i'^i^iT have unveiled many fea- 
tures which stand in stark contrast, e.g., to the Halperin- 
Saslow (hydrodynamic) picture^ of a linear dispersion re- 
lation for low-energy excitations. 

Here, we shall assume that the spin magnitude is |S| » 
1 , so that quantum fiuctuations can be safely neglected^ 



(classical limit). 

An alternative to using the Mattis picture can be pur- 
sued by studying usual spin glasses (i.e. with random ± J 
bonds) in the high-field limit, as this additional feature 
stabilizes a ferromagnetic-Hke ground state while still in- 
corporating quenched (bond) disordes^iifliii. However, 
results thus obtained differ rather drastically from those 
pertaining to the zero-field case. In fact, it has been 
found that, even in zero field and space dimensionality 
d — 1 where frustration effects are absent, "unmagne- 
tized" spin glasses (i.e. in which the concentrations of 
ferro- (p) and antiferromagnetic {1 —p) bonds are equal) 
differ substantially from their "magnetized" (p ^ 1/2) 
counterpart! 



.12 



In this paper, we investigate the spectral and local- 
ization properties of Heisenberg-Mattis spin glasses. Our 
emphasis is on unmagnetized systems in space dimension- 
alities d = 2 and 3, at T = 0. We use numerical transfer- 
matrix methods to calculate Lyapunov exponentsi^iiii^, 
and eigenvalue-counting theorems, coupled with Gaus- 
sian elimination algorithmaiS^ii, to evaluate densities of 
states. Though early numerical studies^ already high- 
lighted a number of distinctive features exhibited by such 
systems, motivation for further research is to be found 
in recent theoretical insights£»iS, especially in connection 
with the low-energy, long-wavelength regime. 

In Section|n]we recall pertinent aspects of Heisenberg- 
Mattis spin glasses. Section ITTTI renorts on an extension, 
to d = 2 and 3, of the analytical scaling techniques intro- 
duced in Ref.lglfor d — 1; in Section HVl we report numer- 
ical calculations of Lyapunov exponents and of densities 
of states, for d = 2 and 3. Finally, in section conclud- 
ing remarks are made. 
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II. HEISENBERG-MATTIS SPIN GLASSES 

We consider Heisenberg spins on sites of a square, or 
simple-cubic, lattice, with nearest-neighbor couplings: 

n = -Y^ J., s. • s, (1) 

The bonds are randomly taken from a quenched, binary 
probability distribution, 

P{.h,) - p 5{J,j - Jo) + (1 - p) 5{J^3 + Jo) : (2) 

so for p = 1/2 one has the unmagnetized spin glass. 

The Mattis model ascribes disorder to sites rather than 
bonds {Jij Jo Ci 0)) so that the Hamiltonian reads: 

"Hm = — Jo ^ Ci ' ^3 ' (3) 

where Ci — +1 (^1) with probability p (1— p). This way, 
the overall energy is minimized by making — QS, 
which constitutes a (classical) ground state of the Hamil- 
tonian Eq. (El, to be referred to as | 0). Thus, disorder is 
effectively removed from static properties, but not from 
the dynamics, because of the handedness of Heisenberg 
spin commutation relations. Indeed, considering low- 
energy excitations, the equations of motion for the spins 
are, with ft = 1: 

I dSr /dt = ^ Jo G - S- 5f ) , (4) 

3 

where Sf^ ^ Sf ±i Sf etc and j are nearest neighbors of 
site i. So, putting Vi = Q , one gets^, upon application 
of Eq. Q to I 0) : 

iC,idui/dt^'^jQ {ui~Uj) . (5) 

3 

where the Ui are Mattis-transformed local (on-site) spin- 
wave amplitudes. For the eigenmodes with frequency lo 
(in units of the exchange constant Jo), Eq. Q leads to 

ujQiUi='^^{ui-Uj) . (6) 

3 

Goldstone modes are expected to occur, since disorder 
does not destroy the symmetry of the system in spin 
spac^. The relationship of frequency to wave number, 
k, at low energies is characterized by the dynamic expo- 
nent z: 

w oc fc^ . (7) 

In 0? = 1, where the scattering length coincides with the 
localization length^, the definition of k is unique. In- 
deed, numerical calculation a^fti^^i^i? of the d = 1 density 
of states and of the Lyapunov exponent point to the same 
value z — 3/2, predicted analytically-. For > 1 this 
degeneracy is expected to be lifted. As we shall see be- 
low, different exponents come up, depending on whether 
localization or density-of-states properties are being con- 
sidered. 



III. SCALING 

We briefly review the treatment of one-dimensional 
systems, given in Ref. 0. In this case, Eq. (jSJ becomes 

(2 - Ci w) Ui = + Ui+i . (8) 

A transfer-matrix (TM) approachi^ii^, can be formu- 
lated, givingSiiS^: 

f2-QuJ -l\ f Ui \ rr. f ^ f Ui \ 

[u.) = [ 1 j U-J ^ '^"H--J ■ 

(9) 

The allowed frequencies for a chain with N spins and pe- 
riodic boundary conditions, mat+i = ui, are determined 
by det (A at — 1) = 0, where 

N 

An{lj) = Y[t,{uj) ; (10) 

i=l 

equivalently, the condition TrAAr = 2 determines the 
eigenfrequencies. Scaling the system by a linear dilation 
factor b, the dynamics is preserved if the frequencies are 
transformed (w uj'), in such a way that 

TrAAr(cc;) =TrAAr/b(c^') • (H) 

Using properties of the matrices Ti{uj), one finds^ that 
the first-order term (in uj) of Tr AAr(ct;) has a coefficient 
equal to N J2i=i d- Therefore, correspondence of the 
{Ci} with an unbiased random- walk makes the determin- 
ing variable luN^^"^, so that the (length) scaHng of the 
frequencies \s w' — uj 6'^/^ , and the low-energy disper- 
sion relation Eq. ^ has an anomalous power (dynamic 
exponent) z — 3/2. In fact, careful consideration of 
higher-order terms^ shows that the combination N^/"^ lo 
is present to all orders, thus scaling is expected to hold 
even away from the w ^ region (though not the single 
power-law form, Eq. Q). 

A suitable framework for extensions of this treatment 
to space dimensionalities c? > 1 is found in quasi- one di- 
mensional geometries, i.e. L'^^^xN systems with TV » 1. 
In what follows, we shall always make use of periodic 
boundary conditions across the d — 1 transverse direc- 
tions. 

Considering d — 2 for simplicity, a TM can be set up 
on a strip of width L sites, so an L- component vector 
Ui — (uii, ■ ■ ■ ,ULi) corresponds to each column i along 
the strip, with the recursion relation 

where 

T?'(^) = i^/ ^a-u;h, (13) 
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/ being the L x L identity matrix, while a and bi are 
given by: 



The result is: 
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Hence, 



(14) 



Generalizations to higher c? are immediate, with the vec- 
tor Ui now having L'^~^ components, and the matrices 
/, a and b, being L'^-'^ x L-^^^. The (2i'^-i x 21'^-'^) 
matrix Ti is symplectic, that is, its eigenvalues occur in 
pairs {vi, i'^^}: i = 1,. ■ ■ ,L'^~^. Note that matrix A is 
symplectic as well. 

For d > 1, a feature which does not occur in the one- 
dimensional case is that there are transverse momentum 
modes. Returning to d = 2 for illustration, these are in- 
deed the eigenmodes of matrix a in Eq. II14II , with corre- 
sponding energies Sp = 4— 2 cos 27rp/L, p — 0,1, ... , L—1. 

We briefly make contact with the analogous case of 
a homogeneous system of length N, for which b — I, 
[a, b] = 0, and the eigenstates of a are also eigenstates 



of the full hamiltonian, with ojp 



2cos2TTq/N, 



q = 0,1, . . . , N — 1. This reminds us that, in d > 1, the 
energy of a mode is not related only to its longitudinal 
wavevector, as is the case in d = 1. Upon introduction of 
randomness, the commutation relation is destroyed (con- 
trary to the one-dimensional case where both a and b 
are numbers) and, consequently, the interplay between 
frequency- and wavevector- aspects can only be measured 
via the accumulated statistics of many local realizations 
of disorder. Therefore, in d > 1 one may expect the 
picture of a single length controlling both (spatial) at- 
tenuation and (time) oscillation damping^, which holds 
for d = 1 spin glasses, to be replaced by one where each 
of these properties is governed by a distinct quantity. 

We now return to spin glasses. Prom the eigenvec- 
tors of a, "spinor" generalizations can be built, which are 
eigenvectors of A, with eigenvalues (z/p, fj7^) indexed by 
p; one can show that I'p + = £p. While such spinors 
are obviously not eigenvectors of B, the contribution 
given by each diagonal element of Tf^^{u)), correspond- 
ing to fixed p, to the trace of A^(ci;) = Hi^i '^i'^i'^)^ can 
be worked out to first order in lu. Use is made of the fact 
that, analogously to the d = 1 case^. 



JV 



llT^^{cj) = A^ 



1^1 



N 

^=1 



^BoA 



N-e 



0{lo') . (16) 



N 



.N , -N 
P 



Tr(p)nTrH-< 

\^ EE7C^^ + ^("')- (17) 

where Tr {p) denotes the joint contribution of both eigen- 
spinors of A indexed by p (associated respectively to 
eigenvalues Vp and Vp^). 

The "critical" (large scale) behavior is associated with 
small p, in which case Vp, 1, and Eq. II17|I turns 

into: 

N N L 

Tr {p) n Tf^iuj) ^2-^7V5]-^C™^ + 0{lo^) . 

i=l 1=1 m=l 

(18) 

One can readily see that, for generic d > 1, this trans- 
lates into: 

AT 

Tr(p)[]T,(u;)^2- 

1=1 

N /d-1 ^ Li \ 
-^NY,{\{- ^ Cm,...m,_A+0{u^) . (19) 
1=1 \i=l ' mi = l ) 

In the second term of Eq. (jl9ll . one has a sum of iV x ii x 
• ■ • X Ld-\ binary random variables, so this is gaussian 
distributed with rms value: 



/ ^ 1 \ 1/2 



nti ^ 



1/2 



(20) 



Upon scaling of linear dimensions by a factor 6, under 
which frequency scales as w — > 6^ w, and requiring invari- 
ance of the term given in Eq. lf2(Hl (see Eq. ifTTll V one 
gets: 



d 
2 




(21) 



Consideration of the terms in Eq. IjlHIl of higher than first 
order in uj shows that, unlike the d = 1 case, the trace 
of the full TM is not just a function of the variable given 
in Eq. ll2?Hl . because complicated sums occur, involving 
both longitudinal and transverse wave vectors. 

This is in line with the reasoning presented above, to 
the effect that the simultaneous presence of both longi- 
tudinal and transverse degrees of freedom invalidates the 
single-length picture, predicted analytically^ and numer- 
ically confirmedifliiSiii for d = 1. 

While it is plausible to expect that, for some low- 
energy regime in d > 1 the scaling result, Eq. lf2T|l might 
hold true, direct verification is called for. 
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Figure 1: Scaling plot for localization lengths on strips of 
a d = 2 system, against to , with z = 1 as predicted in 
Eq. 12111 . Strip widths L as given by symbols. Line corre- 
sponds to 2/ (X x~^^^ and is a guide to the eye, showing how 
the effective d = 1 regime sets in for very low energies. Insert: 
unsealed data for L = 4 and 14. 



Figure 2: Scaling plot for localization lengths on bars of a 
d = 3 system, against uL^, with z — 1/2 as predicted in 
Eq. pH . Bar cross-sections are L x L, with L as given by 
symbols. Line corresponds to j/ cx x~^^^ and is a guide to the 
eye, showing how the effective d = 1 regime sets in for very 
low energies. Insert: unsealed data for L = 4 and 10. 



IV. NUMERICAL ANALYSIS 
A. Lyapunov exponents 

The procedure for calculating Lyapunov exponents on 
strips or bars is the same as that used for Anderson lo- 
caHzation problema^^. Indeed, in both cases the TM is 
symplectic, and one can use Oseledec's theorem and dy- 
namic filtration to extract the smallest Lyapunov expo- 
nent, whose inverse is the largest localization length. For 
Heisenberg spin-glass chains, this has been donei^ii^, nu- 
merically confirming the result z = 3/2 obtained analyt- 
ically in Ref. 0. 

We have investigated strips of widths L — 4,6, ■ ■ ■ ,14 
in d = 2, in which for each energy uj we took TV = 10^ 
iterations of the TM, and bars with L x L cross-section, 
L = 4, 6, 8, 10 in d = 3. In d = 3 we used N = 10^ for 
L = 4, 6, 5 X 10^ for L = 8, and 1 x 10^ for L = 10. 

In contrast with d = 1, here one must take into account 
finite-size effects, introduced via the transverse dimen- 
sion L, thus calculated locaHzation lengths are denoted 
by Al. Using standard finite-size scaling theorjii^, it is 
expected that the behavior of scaled localization lengths 
Xl/L, when plotted against uL^, will allow one to infer 
the bulk {L — » oo) properties of the system. 

In Fig. n] we see that in c? = 2 good data collapse, 
extending as far as x = ui ~ 0.3, is achieved when z = 
1, as predicted in Eq. 1I2TI1 . At the low-energy end, x < 
0.03, the quasi-one dimensional character of the strips 



begins to dominate, and the scaling curve crosses over to 
the effective d = 1 regime characterized by w ~ k^^^ . 

In Fig. 121 the scaling plot for d = 3, with z — 1/2 
as predicted in Eq. (l2T|l . is exhibited. The quality of 
data collapse is remarkably inferior to that of d = 2 
data. An examination of the behavior of Xl/L against uj 
shows that curves corresponding to pairs L,L — 2 have 
well-defined crossings at low energies u < 0.05. The 
usual interpretation of these, in the finite-size scaling con- 
text, would point to a localization-delocaHzation transi- 
tio n^i?!^!? . However, we have found that the locations of 
crossings appear to approach uj = with increasing L. 
This would be consistent with the idea that all magnons 
are delocaHzed in d = 3, which is supported, e.g., by 
the field-theoretical results of Ref. [3. We postpone a dis- 
cussion of this point (and similar ones associated to the 
behavior found above for d = 2) , to Section 

B. Densities of states 

The calculation of densities of states per unit en- 
ergy interval (DOS), V{uj), and their integrated coun- 
terparts (IDOS), N{lj) = J';^^Viuj')duj', makes use of 
eigenvalue-counting theorems^SiSiiS^. Our implementa- 
tion resorts to Gaussian elimination algorithms on quasi 
one-dimensional geometries {L'^^^ x N, with TV ^ L), 
and closely follows the steps described in Refs. 'ifiilTl 
where the systems under investigation were, respectively. 
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Figure 3: Double-logarithmic plot of density of states T>(yi) for 
spin-glass chain, calculated by Gaussian elimination. Chain 
length A*' = 10^ sites. Thick line is the exact Derrida-Gardner 
resuHt^i (with coefficient doubled, on account of different nor- 
malization, see text). 



phonons in disordered solids, and tight-binding electrons 
(Anderson localization) . The key feature shared between 
these problems and the one studied here is the fact that, 
for an L'^~^ x iV system with periodic boundary condi- 
tions across, the hamiltonian has a (2 x L'^~^ -|- 1)- diag- 
onal form, i.e., it can only have non-zero elements in the 
L"^^^ lines above, and L'^~^ Hues below, the diagonal. 

We consider the characteristic matrix, which in the 
present case is C = (oj I — Ti, where (uj I is a diago- 
nal matrix with [(^uj I]jj = (j uj {j = site index), for an 
l^d-i ^ system. Evaluation of its diagonal elements 
via Gaussian elimination enables one to obtain the IDOS 
for any energ y^^'^T , thus the DOS may be calculated by 
numerical differentiation. 

For d = 1 the eigenvalue counting (used., e.g., in 
Ref. IT^ , may, alternatively, proceed via enumeration of 
nodes of the amplitude ratios which enter the evalua- 
tion of the (single) Lyapunov exponen^iSii^. In order to 
test our recursion and elimination algorithms, we applied 
them to this case and compared the outcome with that 
from node-enumeration . Results are identical to within 
numerical accuracy, and the set produced by Gaussian 
elimination is depicted in Fig. By sampling ener- 
gies separated by logarithmically uniform intervals, we 
achieved a detailed view of the uj ^ region, which is 
difficult to isolate in the corresponding DOS results of 
Refs. Hclll^lT^ (where linear binning was used). One sees 
that the relationship I?(w) oc u!~^^^ is valid for more 
than two orders of magnitude in energy, up to w ~ 0.3. 
For guidance, we have also included the exact Derrida- 



Gardner results^. Since we have considered only positive 
energy states in our calculation, the appropri ate p ropor- 
tionality coefficient is twice that given in Ref. I23. 

For higher-dimensional cases, is is worth mentioning 
that the algorithms used here are much less computa- 
tionally intensive than their Lyapunov-exponent counter- 
parts. For an L"*-^ x TV system, the computational time 
rises as L'^^'^~^^ x N for the formepiS, and approximately 
as L^'"*^^^ X N for the latter. This is mainly because 
of the frequent mutual orthogonalization of 2 x L'^~^ 
iterated vectors, which is necessary in order to avoid 
cross-contamination between eigenvectors associated to 
different Lyapunov exponents. Therefore, for DOS and 
IDOS it is usually possible (except for very low energies in 
d = 3, see below) to work with systems whose transverse 
dimensions L are large enough that finite-size effects are 
of little import. It remains only to make sure that the 
sample length N is long enough, in order to achieve ad- 
equate sampHng of quenched disorder configurations. 

We examined the effect of finite transverse dimensions, 
by evaluating pure-system quantities and comparing our 
results to the exact ones. Though, having zero net mag- 
netization, the spin glasses studied here are closer to an- 
tiferromagnets (AF) than to homogeneous ferromagnets 
(FM), the DOS and IDOS of magnons in the latter ex- 
hibit some distinctive features, whose numerical repro- 
duction is a non-trivial test of the adequacy and accuracy 
of our methods. For FM in d = 2, already with L — 25, 
N = 2500 the IDOS is at most 3% off the exact value. 
This largest discrepancy happens close to a; = 4 where 
the analytical IDOS exhibits an infiection point, on ac- 
count of the DOS's logarithmic Van Hove singularity at 
the band center. Increasing L or N does not significantly 
reduce the deviation close to uj — 4; however, it does 
improve agreement elsewhere on the energy axis. The 
calculated DOS is rather sensitive to discrete-lattice ef- 
fects; nevertheless, the consequent oscillations are again 
much diminished by increasing L, N. For d — 3 FM, the 
relatively featureless IDOS is easier to reproduce. With 
L = 16, iV = 25600, deviations are down to, at most, 
1.5% (though the DOS still displays somewhat large os- 
cillations, especially around the "knees" aX uj — A and 8). 
Fig. 01 shows representative results, which are useful as 
guidelines for the investigation of disordered systems in 
d — 2 and 3 via Gaussian ehmination. 

Turning to pure AF systems, for which the respective 
bandwidths are wff^ =4 (d = 2), and 6 {d = 3), again 
relatively small transverse dimensions L provide results 
which closely follow the analytic values, except at very 
low u!. In this limit, the fact that the finite L quan- 
tizes the transverse momentum leads to effective one- 
dimensional behavior {T>{uj) ^ lu^, N{uj) ^ uj^) for uj 
less than a crossover frequency uJm = AAF{d)/L^, z = 1. 
With the units used in this work, we found Aaf(2) — 12, 
v4af(3) — 20. The effect is more pronounced here than 
for FM, where z = 2 and, consequently, the onset of this 
sort of behavior occurs at much lower energies. Fig. El 
highlights the worst case of d = 3. For completeness, the 
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Figure 4: DOS {V{lu)) and IDOS {N{uj)) for pure FM sys- 
tems against normalized energy cj/ojq: analytical (lines) and 
calculated by Gaussian elimination (points). Dashed lines 
and triangles: ^{u); full lines and squares: N{lo). (a): d — 2, 
L = 100, N = 50000. (b): d = 3, L = 16, = 25600. 
The pure-system (FM) bandwidth is ujo —8 [d — 2), and 12 
(d = 3). 



Figure 5: Low-energy IDOS, N{uj), for pure AF in d = 3 
against energy uj: analytical (full line) and calculated by 
Gaussian elimination on x N systems, A'^ = 500 (points, 
connected by dashed lines). Triangles: L = 16; squares: 
L = 20; circles: L = 24. Inset: full-band IDOS (same axes as 
main Figure). Analytical (full line), and Gaussian elimination 
with L = 16 (triangles). 



inset of Fig. El shows that, even for L = 16 where these 
low-energy discrepancies are rather severe, agreement 
with analytical forms is quite satisfactory elsewhere. 

We now return to disordered systems. In Fig.El results 
for the Mattis spin glass in d = 2 are presented. We used 
L = 250, = 2.5 X 10^. The number of sites entering the 
calculation was more than one order of magnitude larger 
than in that for a pure FM, whose result is exhibited 
in Fig. 01 (a). From examination of shorter runs for the 
disordered case, it appears that the features displayed in 
Fig.Elare rather stable and well-converged. For this value 
of L, the crossover to one-dimensional behavior, referred 
to above, is confined to w < 0.05, leaving a broad win- 
dow at low energies for which genuine two-dimensional 
behavior can be observed. The main distinctions of the 
IDOS from its pure-system (FM and AF) counterparts 
are: (i) close to ui — 4, the upper limit of the AF band, 
the FM IDOS's inflection point is replaced by a seem- 
ing "knee", with a short flat section; and (ii) saturation 
is reached below the FM band edge uiq = 8, but above 
the AF edge uj^^ = 4; by w = 6.7 the IDOS is aheady 
within less than 1% of unity. Similar effects can be seen 
in early numerical work^, though in that Reference sat- 
uration appears to be reached only above the FM band 
edge, at a; ~ 9.0. 

It has been predicted^ that, since d = 2 is the crit- 
ical dimensionality in this casein, the two-dimensional 
spin glass will behave as a pure (AF) system (namely. 
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Figure 6: DOS CDico)) and IDOS {N{uj)) for Mattis spin 
glass in d = 2, against energy uj, calculated by Gaussian 
elimination. Triangles: T>{lu); squares: N{uj). L = 250, 
iV = 2.5 X 10^ 
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V{uj) - w\ N{uj) - u)'^), with logarithmic corrections. 
At low frequencies, the real part of the dispersion rela- 
tion is expected to follow the expression^ii: 

Re w oc ^ , (22) 

where A is a momentum cutoff, reciprocal to the mini- 
mum wavelength of magnons. From Eq. II22|I . one can 
work out the predicted behavior of the IDOS at low en- 
ergies. This turns out to be: 

7V(w) oc In , (23) 

where is a cutoff frequency, corresponding to the mo- 
mentum cutoff A. 

We have tested the prediction, Eq. II23II . against our 
data, with the results shown in Fig. [3 A fit of the raw 
data (crosses in Fig.CJ to pure power-law behavior gives 
N{lu) ~ w^, with the effective exponent x ~ 1.62. On the 
other hand, plotting N {to) / hi{fl / uj) against (squares 
in Fig.CJ removes just about all the curvature, provided 
that a suitable value of ft is used. A linear least-squares 
fit of data for 0.05 < uj < 0.5 (shown as a full line in 
Fig. gives fl — 5.8(1), broadly consistent with the ef- 
fective bandwidth > 6.7 found above . Keeping SI — 5.8, 
and fitting N{uj)/ Inlfl/uj) to a power law dependence 
over the full interval 0.05 < uj < 1.0, would give an effec- 
tive power X ~ 1.04. 

We undertook similar calculations for the Mattis spin 
glass in d = 3. Since one is above the critical dimension- 
ality in this caseSiifi, the three-dimensional spin glass is 
expected to behave as a pure (AF) system, at least at 
low energies and long wavelengths (namely, ^{lu) uj'^, 

Similarly to the pure d — 3 AF, for the ranges of 
L within relatively easy reach of our calculations, the 
low-frequency spectrum exhibits a crossover towards one- 
dimensional behavior. With the terminology introduced 
above, this happens for u < Um, = Ascid)/ L; by 
examining the sequence L = 16, 24, 30, 36, we estimate 
Asg(3) — 11, just over half the corresponding value for 
pure AF. Thus, such effects are once more confined to low 
energies. We have found that, for uj > 1.2, the L — 16 
curve is within less than 3% of those corresponding to 
larger L, which are grouped together even more tightly. 
Fig. IHl presents an overall picture of results, for L = 16, 
N = 2.56 X 10^. Again, early saturation occurs. The 
IDOS is within 0.1% of unity by w = 9.4, just over three- 
quarters of the FM band width ujq — 12. A kink, similar 
to the one occurring in d = 2 but less intense, arises close 
to the center of the FM band (and top of the AF one), 
Lu — 6. Both features show up in Ref. 0, though with 
saturation occurring at a slightly higher energy (but still 
within the FM band). 

The low-energy behavior is shown in Fig.|2| For L = 36 
we have found that least-squares fits of our calculated 




Figure 7: IDOS {N{uj)) for Mattis spin glass in d = 2, for low 
energies, against to^, calculated by Gaussian elimination. L = 
250, iV = 2.5 X 10^ Crosses: N{lu); squares: N(iu)/ \n{n/uj), 
with — 5.S [see Eq. ^2',^ . and text]. Full line is a linear 
least-squares fit to data for 0.05 < uj < 0.5. 




Figure 8: DOS {V{uj)) and IDOS {N{uj)) for Mattis spin 
glass in d = 3, against energy uj, calculated by Gaussian 
elimination. Triangles: T>{uj); squares: N{uj). L — 16, 
iV = 2.56 X 10^ 
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Figure 9: Double-logarithmic plot of IDOS (A^(t^)) for Mattis 
spin glass in d = 3, for low energies, against energy cj, calcu- 
lated by Gaussian elimination. Squares: calculated points, 
L = 36, iV = 6.48 X lO'^. The straight line is a power- 
law fit with slope 2.97 (from least-squares fit of data for 
0.45 < 1^ < 1.0). Also shown is IDOS for pure AF, calcu- 
lated for a cube with 100^ sites. 



data (excluding the very low-energy intervals where one- 
dimensional behavior takes over) give z = 2.97(5), if 
we keep to w < 1.0; including higher energies (e.g. 
uj < 2.5 — 3.0) results in a slight decrease of effective 
exponents, down to z ~ 2.75. On the other hand, fits of 
the numerically-evaluated analytic IDOS for a cube with 
100'^ sites (shown in Fig.EJ, when restricted to w < 1.5, 
give an effective z — 2.82(1); it is only when the upper 
limit is raised to w = 3.0 that one reaches z — 2.99(1). 
This is because, in the low-energy Hmit, discrete-lattice 
effects still persist, which induce slight deviations of ef- 
fective behavior away from the exact value z = 3. In 
summary, it is only in the very low-energy limit a; < 1.0 
that the d — i SG N(ui) indeed exhibits the lo^ depen- 
dence characteristic of the pure AF. 

Therefore, we conclude that our low-energy data 
are consistent with the indications of Refs. »71il^ that 
magnons in the d — "i Mattis SG display the same low- 
energy behavior as in a pure AF. However, the respec- 
tive amplitudes differ, as is apparent by the roughly con- 
stant distance between SG and AF data in Fig. |2l Writ- 
ing Nx{lo) = ay^iu^ (X=SG, AF), we get from our fits: 
asc/oAF = 4.1(1). 

A calculation of the amplitudes, along the lines of 
Ref. H yields asc/oAF = [2lf/^ = 5.281..., where 
/ = 1.516386 ... is Watson's integral^i. We believe the 
order-of-magnitude agreement found between our numer- 
ical estimate and this result is satisfactory, given that 



disorder is treated only approximately in the latter ap- 
proach. 

V. DISCUSSION AND CONCLUSIONS 

The preceding results are consistent with our state- 
ment, made in SeclUl that the single-length picture which 
prevails in d = 1 cannot be ported to higher space di- 
mensionalities. In order to make contact with the one- 
dimensional case, we will refer to the indices emerging 
from the analytical scaling of Sec. lIIII and from the Lya- 
punov exponent calculations of Sec. IIV Al as zl, while 
those originating from the results of Sec. IIV Bl (plus the 
relationship N(lo) ~ w'^/^) will be denoted by z^. 

The analytical scaling predictions z^ = 1 (d — 2) , zl — 
1/2 {d = 3), are confirmed by our Lyapunov exponent 
calculations, though the width of the energy intervals for 
which scaling holds is larger for the former (5 x 10~'^ < 
i^L < 0.3) than for the latter (5 x lO^^ < u}L^/^ < 
5 X 10-2). 

In d = 2, the curves of Al/L against uj are essentially 
parallel for uj < 0.1, down to the lowest energies inves- 
tigated; for fixed w, A^/L decreases with increasing L. 
This indicates the absence of a delocalization transition, 
i.e. all modes are localized in d = 2, in agreement with 
Refs. On the other hand, our result zl = 1 implies 

that the localization length diverges at low energies as 
£ioc '-^ uo~^ . This is in contrast with the field-theoretical 
prediction of Ref. 0, according to which ^loc ~ w"^/^^'^. 

For d = 3, as mentioned above, the curves of Xl/L 
against u cross each other at low energies. For the 
(L, L — 2) — (6, 4) pair, the crossing occurs at ~ 0.04, 
while for (10, 8) it moves to lower energy lo ~ 0.015. We 
interpret this as a residual finite-size effect, which will 
properly vanish with increasing L, and see no reason why 
the estabHshed ide^J^ that all excitations are delocal- 
ized in d = 3 should be challenged on the basis of such 
result. 

A connection of our predictions for zl with the litera- 
ture can be made as follows. The analysis of Refs. i^jJi was 
carried out by assuming a well-defined (real) wavevector, 
thus implying the complex dispersion relation: 

uj{k) =ujR{k) +iV{k) . (24) 

On the other hand, our TM formulation gives a speci- 
fied (spatial) amplitude decay ratio A~^ for a fixed (real) 
frequency, which then envisages a complex wavevector, 

k = kR + iki , Xr^kJ^ . (25) 

One can then plug Eq. lf25|l back into Eq. (jjjl, taking 
into account the specific dependencies of ujji and T on fc, 
and force to to be real in the latter. 

For d = 3, one expect8iia;fl(fc) ~ fc, r(fc) ~ A:^, consis- 
tent with small fine broadening at low k (i.e. propagating 
modes). From this, one then gets: 

A-i ~ (d = 3) , (26) 
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so that the scaling variable is indeed ljL^/'^. 

For d = 2, a similar argument can be made (now on 
somewhat flimsier grounds, because all modes are ex- 
pected to be localized, so the real and imaginary parts of 
the dispersion relation may be of the same order of mag- 
nitude). Ignoring logarithmic corrections, the results of 
Refs.00are: ujR{k) ^ fc, r(fc) ~ fc, from which we get: 

\-^^uj (d=2), (27) 

again consistent with the d = 2 scaling variable being 

LOL. 

The outcome of our density-of-states calculations for 
d = 2 can be very closely fitted, for low energies 0.05 < 
u! < 0.5, to the logarithmically-corrected form predicted 
in HvL 7 (see Eqs. (|22l, JSJ, and Fig. 01. Furthermore, 
one gets z^j = 1 plus enhancing logarithmic corrections 
(recall the effective exponent ~ 1.62 from Fig.[Z||, which 
is in line with the vanishing of group velocity (mode soft- 
ening}^ as a; ^ 0. 

Finally, our d — 3 density-of-states results are again 



consistent with the pure AF behavior predicted^iiS to 
hold above dc = 2. Thus we have z^^ = 1 in this case. 
However, the amplitudes of the low-energy power-law be- 
havior differ, and we have found asc/oAF = 4.1(1). 
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